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In a superconductor with magnetic impurities, Kondo scattering results in the formation of local- 
ized states inside the superconducting gap. We show that inelastic electronic transitions involving 
quasiparticle scattering into and out of the localized states may result in significant changes in the 
non equilibrium properties of the superconductor. Using the model of MuUer-Hartmann and Zittartz 
for the extreme dilute limit, and including both deformation potential and spin-lattice coupling we 
have calculated the rates of such inelastic transitions between continuum and discrete states, and 
shown that they may greatly modify quasiparticle interactions. The individual processes are: quasi- 
particle trapping into discrete states, enhanced recombination with localized quasiparticles, and 
pair breaking and detrapping of localized quasiparticles by sub-gap phonons. We find that all these 
processes give rise to clearly distinguishable temperature dependences of the kinetic parameters. 



I. INTRODUCTION 
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• ' The study of the effects of magnetic impurities in superconductors originated with the pioneering work by Abrikosov 

j^ . and Gor'kov^. Recent interest in the subject has been greatly stimulated by direct observation of the states bound to 

C^ impurity atoms^ which has led to extensive experimental and theoretical work in both conventional and unconventional 

I , superconductors^. An important consequence of such intra gap bound states which has not previously been considered 

''O ' is their role in providing enhanced trapping and recombination at impurity atoms in analogy with deep levels in 

G , semiconductors^. Thus, quasiparticles initially in continuum states may undergo inelastic scattering with phonon 

emission and become localized in the vicinity of impurity atoms, which will act as recombination centers and provide 

rapid thermalization of a non equilibrium initial distribution. The formation of an intra gap band of impurity levels, 

possibly even overlapping the ground state, will modify the temperature dependence of thermalization. Finally, 

activation of localized quasiparticles into the continuum spectrum results in an anomalous temperature dependence 

of the observable parameters characterizing the non equilibrium state, such as quasiparticle lifetime"'. Previously 

neither mechanisms of coupling nor transitions between the continuum and discrete states bound to impurities were 

discussed. In this paper we will show that electronic transitions between the continuum and bound states occur both 

due to deformation potential and spin-lattice interaction. As will be described later there is strong evidence that such 

a scenario has already been observed experimentally. 

We will consider the dilute impurity limit c <SC 1. Here c is the dimensionless impurity density in units of the 
QQ ' condensate density 2A^(0)A, where A^(0) is the density of states at the Fermi level per spin in the normal state, and 
f^ . A is the gap. The problem can be explicitly solved within the model originally developed by Miiller-Hartmann and 
Zittartz^ for quantum spins in fully gapped superconductors. In this model the Hamiltonian of the system is taken 
in the form 



H' H = Ho + H' (1) 

where Hq is the Hamiltonian of an ideal superconductor and H' describes the interaction between impurity atoms 
and conduction electrons. The corresponding interaction potential has the fornti 

w(r)=^[iii(r-R,)+M2(r-R«)<T-5] (2) 

i 

where Ri is the coordinate of the impurity atom S is its spin and cr'^g'^ are the spin Pauli matrices. Here the first term 
describes the spin independent part of the impurity scattering potential, and the second term the exchange interaction. 
To consider phonon assisted electronic transitions involving the bound states we must include also terms, describing 
the electron-phonon interaction through both the deformation potential and the spin-lattice interaction; these are 
derived by expanding the first and the second terms in expression (U) respectively to include the displacement of the 
impurity atom from its equilibrium site. In the four-dimensional matrix formalism the full interaction Hamiltonian 



describing phonon assisted electronic transitions has the form 

H,„t^^fdrt^+{r)V{r)^P^v) (3) 

where 4''^{r) and 4'i^) ^-I'o four-component operators and ^(r) is the 4x4 matrix of the form „ »tr I? where 

v{r) — ^j QiV (uiCTo + U2(J2 ■ S), w*'' is the transposed matrix, and Qi is the lattice displacement of the impurity due 
to vibrations. 

In a superconductor described by the Hamiltonian ([T]) for an impurity with antiferromagnetic exchange, bound 
states split off from the gap are formed^. For the dilute limit the shift of the gap edge remains small, being pro- 
portional to the density of impurities. Therefore we may disregard all effects of modifications to the continuum 
quasiparticle spectrum relating to level shift or to broadening. Our objective is to determine the spatially averaged 
Green function for the continuum spectrum, with the interaction given by ([3]). In the limit c <C 1 spatial averaging can 
be carried separately for all elements of the Dyson equation. Indeed, the interaction loop itself will give a contribution 
proportional to the number of discrete states and hence to c. Therefore, replacing the external Green functions by 
spatially averaged ones introduces inaccuracy in terms only of the order of c^ <C 1, so that the only function remaining 
to be averaged is that inside the loop. After spatial averaging, the Fourier transform to the momentum space, and 
analytical continuation from the imaginary to the real axis iujn —^uj + i6, we may write an expression for the self 
energy 



Im 



( G(p',z')l5^,(P,P')P'^o J^(p',^')l.9q,(P,P')P*^2 
1^ -F+(p',z')l.9q,,(p,p')P*^2 -G(p',z')l4,,(P,P')P^o 



tanhl -^ I -|-coth(2|^ 

z — z' -\- i5 



(4) 



Here G(p', z'), G{p\ z'), F{p' , z') and i^+(p', z') are spatially averaged electronic Green functions obtained within the 
model of Miiller-Hartmann and Zittartz in the dilute limit. The coupling strength gq ,(p,p') is introduced through 



V^|2 _ |„(1) (-r. r,'M2 _|_ Q^Q j_ 1M„(2) 1^ „M|2 

h 

We obtain for the rate of quasiparticle transitions from the state (p, e) the following expression 



l4,,(p,p')P = l5i'i(P,P')P±^(^+l)l5S(P'P')l' 

'^q'-^P' p')'' = 2MNu . ^p " p'' ^q^^-)'i"i^2(p - p')P (5) 



rp(e) = ~ 2eZ (O) ^^^P ^ e)ImSi,pft + (Cp - e)ImEi^p,i 



Aim I]2,j,,. + S+ ] (6) 



^2,ph 

where 2'i(0) is the real part of the renormalization parameter. Introducing coupling constants analogous to the 
Eliashberg constant 

d'P f d^P' Y- U(1.2)^„ „/^|2s 



f d^p 

we obtain 



a^i,2) {^)F{z) - " " ^-^p (7) 

is 



^^^y- z 



^J £ dzdz'F{z){al{z)Re 



Gra{z') ' jFrniz') 



+S{S+l)al{z)Re 
{ 



tanh I — 



tanh — - — coth 
2T 



G„,(z') + -F„(z') } 
S{e-z-z') 
S{e + z-z')} 



coth (—] 
\2TJ 



(- 
\2T 



(8) 



Here we have introduced the new notations Gm{z) and Fm{z) for the spatiaUy averaged Green functions. The exact 
expressions for the Gm{z) and Fm{z) have the form 



G™(e) 



e(e) 



Fraie) = 



A(e) 



The energy e and the order parameter A satisfy the following equations^ 

? = e + AEi(g,A); A = A + AE2(e,A) 
where 



(9) 



(10) 



Si(y,A). 
S2(y,A) 



c \/ y — 1 

— 2 _ „2 y^y ~ yo) 



^■^ y - yg 
c \/y^ - 1 
«7r y^ - vl 



yo(y-yo) 



(11) 



Here y = e/A, yo = eo/A < 1 and eo is the discrete intra gap level. To find solutions to the main terms, we may use 
the simplified equations obtained from (|10p by taking the renormalised parameters y = e/A and yo = eo/A in the 
denominators in (jlip . The simplified equations then become 



z — c- 



1 



1 



V_ 

\v\ \y-yo y + yo 
\y\ \y~ya y + yo 



(12) 



where c = c— . In order to solve the simplified equations we first note that from (jlip E2 — —yo/y^i- 

27ryo 

Above the gap edge both jSij ~ c and IE2I ~ c; as pointed out above we ignore these corrections. Inside the gap the 

solution for the imaginary part S" has the form 



E'/ 



c-\{\y\~yoy&(c-^i\y\^yof 



(13) 



Within the range E" 7^ 0, the real part of the self-energy 'E,[ is given by 



si = ^feo-M) + ^ 



(14) 



Outside this range, but still inside the gap, T,'-^ remains finite with a dependence on impurity concentration changing 
from y/c at the edge of the S" ^ range to c away from it. 
Finally, inside the gap we obtain 



i^^G™(y) = (/_^J)°3/2 V^"-i(M-^o)^Q (^- \i\y\-yor 



^-p-iy) = f^^f"- \i\y\-yoy^ (c- liM^yof 



(15) 



These expressions describe the normalised density of bound states inside the gap in a superconductor with magnetic 
impurities. This distribution is sharp, with both its width and height being proportional to ^/c. It is easy to confirm 



that / dyKeGmiv) = c, corresponding to one quasiparticle bound state for each impurity atom. Inside the gap, 

from P^ . we also obtain ReFm{y) — yoReGm(y). 

Using the Green functions given by (jlSp and ^ we may now analyze the different inelastic transitions. Firstly, 
a quasiparticle initially in the continuum state may become trapped. For the trapping rates we obtain 



2T{uj) 



de'(e - e') 



ReG„(e') 



[1 -/(£')] = 



Ttrapie-) Tt- Jo A2 


1 lA A , „ / 1 1 \' 

^ + ^ - _ReF,„(e') ^ - ^ 
Ti T2J e \Ti T2J 


He - e') + 1] 


fi Ae f2 Ae 


He - eo) 


+ i][i-/M] 



(16) 



where n{e) and /(e) are the phonon and quasiparticle distribution functions. The characteristic relaxation times for 
phonon assisted scattering on magnetic impurity in the host lattice tm 2) can be written in the form 



1 1 af(A) /A 



fi roa2(A) yT, 



l_ _ S{S+l)al{/\) /A 
f2 " To a2(A) \t. 



(17) 



where a is the parameter entering the Eliashberg constant, tq is the superconductor characteristic relaxation time 
for deformation potential coupling, and Tc is the critical temperature. The top bars in this notation emphasize that 
these characteristic times are for phonon assisted scattering on a magnetic impurity. An order of magnitude estimate 
of the ratio ro/fi_2 can be obtained by direct evaluation of factors a\ 2- Thus 



To 



UeiJ Aas T2 \UeiJ Aa^ 



(18) 



where ui, U2 and Uei are the characteristic values of electron-impurity, exchange interaction and electron-ion potentials 
respectively, v^ is sound velocity and a^ is the radius of the bound state. 

The recombination rate via a bound state calculated from ([5]) and P^ is given by 



ReG„,(e')(^ + ^ 

n T2 



Ke + e') + l]/(e') = c 



rfl,t(e) 



de'{e + e') 
A2 



-ReF,„(6') {^-^ 

e \Tl T2 



1 {e + eof 



1 e2-e2l 



Ti Ae r2 Ae 

[n(e + eo) + l]/(eo) 



(19) 



The expression can be written in a more familiar form by introducing the appropriate recombination coefficient i?f 
and density of trapped quasiparticles rit 



T ji,t{e) ^ Rtnt] Rt 



1 



2iV(0)A 



1 {e + e^f 



1 e' 



Tl 



Ae 



7=2 Ae 



(20) 



which describes the maximum recombination rate in the absence of a phonon bottle-neck effect. Comparing the 
recombination coefficient on traps Rf with that in ideal superconductor R we obtain 



R 



Uel / V Uel 



Atts 



(21) 



Taking values for (ui/uei) , (^2/^61) — 1, Vs— 3T0^cm/s, A = 0.5meV, S=l and Os^lnm we obtain Rt/R >10. 
This ratio indicates the dominance of recombination at impurities due to the larger magnitude of the spin-lattice and 
deformation potential coupling constant with discrete levels originating from the appearance in ([5]) of a form-factor for 
phonon assisted impurity scattering, instead of the momentum conservation law. Similarly, comparing the maximum 



recombination rates under quasi equilibrium conditions for the two different processes, we obtain 



RtntT 




(22) 



where rit^T and tit are thermal distributions of trapped and mobile quasiparticles. Hence, even for a small impurity 
density, recombination on the traps at low temperatures is a stronger process because of the presence of the exponential 
factor. The presence of this factor significantly accelerates recombination at low temperatures in superconductors 
containing concentrations of magnetic impurities which are below trace levels. Moreover, the possible formation of an 
intra gap band of bound states, and also of discrete bound states in the vicinity of the Fermi level, can significantly 
change the observed temperature dependence of recombination and thermalization rates. In some situations the rates 
in an impure superconductor may remain finite even at T = 0. 

In spite of the stronger coupling constant the pair breaking by sub-gap phonons at small impurity density may 
be less efficient than for transitions in the continuum spectrum. In the latter case strong pair-breaking is known to 
slow down the recombination rate because of phonon bottle-necking. Thus less efficient phonon bottle-necking also 
enhances recombination at impurities. The pair breaking rate can be calculated using the Green functions calculated 
earlier to give 



Tpb{n) 



2ro 



n-A 



""Tp/i Jo 

1 

+ 
n 



^{ (^ + ^\ RcG„,(e)ReG,„(r! - e) 
^\ ReF,„(e)Re^™(17 - e)} = c—ri{^) 

T2 J Tph 



(23) 



where Tph is the characteristic pair breaking time of an ideal superconductor and 






Vx(f^:OT)+2\/c 



^ + (1-2^)^ 



(l + j/o)^/^ 
^1 + n/A - yo 



Vx(f2,ao)-2V5e(x(n.yo)-2Ve) 



dWAc-ixin,yo)-t^y 



(24) 



where x(r2, j/o) = f^/A — 1 — j/o and O is the step function. 

The rate of de-trapping from the localized state can be calculated without the need for direct evaluation of the 
broadening of the bound state due to spin-lattice interaction. Since we are interested in the de-trapping rate due to 
transitions into all available states, we may balance scattering into and out of the bound states at thermal equilibrium. 
The result is 



1 



'^de—trap 



1 

71 



zexp (— z^) 



1 



1 



H^ 



(25) 



The formulas ([16]) , p9)) , (|23|) , and (p5|) may be compared with the well known results for an ideal superconductor—. 

Experimental data indicating the possible presence of such processes has appeared previously in several works. 
Thermalization which is several orders of magnitude faster than expected is routinely seen in a number of super- 
conducting absorber materials used in low temperature single photon detectors. In our own earlier experiments on 
non stationary, non equilibrium quasiparticle distributions in superconducting tunnel junctions (STJs) we found that 
many detailed features of the experiments could only be successfully modelled on the assumption of local trapping 
states with well defined energy levels^. However, in the absence of a microscopic model for the trap, the approach was 
purely phenomenological. All inelastic electronic transitions into and out of the traps were modelled with parameters 
determined from fitting to experimental dataS. In experiments on nominally pure Nb, Ta and Al STJs and various 
proximiscd bi-laycr structures, we showed that the observed behaviour was consistent with the local traps being either 
macroscopic regions of suppressed gap as previously seen in low temperature SEM scans^, or due to the presence of 
magnetic impurities. 

Recently the anomalous temperature dependence of quasiparticle lifetimes in Ta and Al films similar to that 
observed in^ was reported by R.Barends et. al detected through measurements of kinetic inductance and the obser- 
vation of noise spectra^*^. An important feature of the results obtained by this technique was the spatial homogeneity 



of the response, indicating the intrinsic character of the traps and suggesting that they were associated with magnetic 
impurities distributed homogeneously throughout the superconductor. The new mechanisms for inelastic scattering 
of the electrons may also be important for electron decoherence in normal metals with Kondo impurities, currently 
the subject of great interestiiii^ . 

We acknowledge valuable discussions with T.M.Klapwijk, R.Barends and J.R.Gao. 



A.A.Abrikosov and L.P.Gor'kov, Sov.Phys. JETP, 12, 1243 (1961) 

A.Yazdani, B.A.Jones, C.P.Lutz, M.F.Crommie, and D.M.Eigler, Science, 275, 1767 (1997) 

A.V.Balatsky, I.Vechter and Jian-Xin Zhy, Reviews of modern physics, 78, 373 (2006) 

A.Poelaert, A.G.Kozorezov, A. Peacock, J.K.Wigmore, PRL, 82,1257 (1999) 
^ A.G.Kozorezov, J.K.Wigmore, A. Peacock, A.Poelaert, P.Verhoeve, R.den Hartog, G.Brammertz, Applied Physics Letters, 

78, 3654 (2001) 
^ E.Miiller-Hartmann and J.Zittartz, Phys. Rev. Lett. 26, 428 (1971) 

^ S.B.Kaplan, C.C. Chi, D.N.Langenberg , J.J.Chang , S.Yafarey, and D.Scalapino, Phys. Rev.Bll, 4854 (1976) 
® A.G.Kozorezov, R.A. Hijmering, G.Brammertz, J.K.Wigmore, A. Peacock, D.Martin, P.Verhoeve, A. A. Golubov and H. 

Rogalla, Phys. Rev. B 77, 014501 (2008) 
^ R.Gross, D. Koelle, Rep. Prog. Phys. 57, 651 (1994) 

° R.Barends, J.J.A.Baselmans, S.J. C.Yates, J.R.Gao, J.N.Hovenier, and T.M.Klapwijk, submitted to PRL 
^ Gassem M. Alzoubi and Norman O. Birge, PRL 97, 226803 (2006) 
^ F. Mallet, J. Ericsson, D. Mailly, S. Unliibayir, D. Renter, A. Melnikov, A. D. Wieck, T. Micklitz, A. Rosch, T. A. Costi, L. 

Saminadayar, and C. Bauerle, PRL, 97, 226804 (2006) 



